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ABSTRACT 

We study a stochastic lattice gas of particles undergoing asymmetric diffusion 
in two dimensions. Transitions between a low-density uniform phase and high- 
density non-uniform phases characterized by localized or extended structure are 
found. We develop a mean-field theory which relates coarse-grained parameters to 
microscopic ones. Detailed predictions for finite-size (L) scaling and density profiles 
agree excellently with simulations. Unusual large-L behavior of the transition point 
parallel to that of self-organized sandpile models is found. 
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There is much interest in the studies of lattice gas systems evolving under conser- 
vative stochastic dynamical rules [0]. When the rules do not satisfy detailed balance 
condition, intriguing properties atypical of thermal equilibrium arise Q. Their map- 
pings onto surface growth models and directed polymer in a random medium|| 
generates further interests and development. Very recently, such connections have 
been exploited to shed light on the unbinding transitions of directed polymer jjj. 

One particularly simple class of lattice gas models consists of particles under- 
going biased diffusion with no interaction except hard-core exclusion. Generally 
known as the asymmetric simple exclusion process (ASEP)|U, it has attracted the 
interests of both physicists and mathematicians for its novel phase transitions @ and 
shock structures Q under various boundary conditions. Generalizing the ASEP to 
a model having two species biased along opposite directions in an L x x L y system, 
Schmittmann et al.^Q observed a transition between a uniform and a non-uniform 
phase as the number of particles is varied, with the latter phase characterized by 
a compact strip of aggregated particles. However, the fundamental question of 
the existence of the transition in the thermodynamic limit has not been answered. 
This is not a simple issue to settle, since the limiting process depends crucially on 
anisotropic finite-size effects induced by making one direction special. The precise 
form of anisotropies is unknown. Nevertheless, a hint is provided by analogous 
continuous phase transitions, which typically have ensemble averages depending on 
L^/Ly with p ^ Depending on how L x — > oo and L y — > oo are taken, undesired 
singularities may arise from L x /L y — > oo or 0. If p is not known, it is not possible 
to proceed. In another context, a cellular automaton model with orthogonal biases 
has been formulated to mimic cross traffic flow||. A similar jamming transition 
has also been found. Whether the transition survives in the thermodynamic limit 
is likewise not investigated. 

Motivated by these open questions, we study in this letter a three-state model 
with two inter-penetrating ASEP orthogonal to each other. By biasing the first 
(second) species along +y (+x) direction, isotropic finite-size effects are ensured (i.e., 
effectively p = 1 above) so that no ambiguity is associated with the thermodynamic 
limit, whereas the essential simplicity of previous models is maintained. Specifically, 
we consider a model on a two-dimensional (2D) square lattice of square geometry 
L x L, with periodic boundary conditions. Occupation numbers {n X;V } specify the 
configuration, where n x ^ y £ {0,1,2} represents a vacancy, a type-1 or a type-2 
particle at site (x,y), respectively. Except hard-core exclusion, the particles are 
non-interacting. In simulations, a site and one of its nearest neighbors are randomly 
picked in turn. If and only if one of them is vacant, a jump occurs with rate: p for 
a type-1 particle hopping along +y, or a type-2 particle hopping along +x; q for 
type-1 along —y, or type-2 along — x; and r for type-1 along ±x, or type-2 along ±y. 
L 2 such trials constitute one unit of time (one sweep). Note that choosing p > q > 
introduces the asymmetry, whereas r > provides transverse diffusions necessary 
for ergodicity (contrary to the traffic flow automaton^ which has r = 0). Using 
Metropolis rules and a parametrization by an "electric field" £ 0, the rates become 
p = 1/4 = r , q = e~ £ /4. 

Simulation results: Simulations are done on wide ranges of parameters: 32 < L < 
256, 0.2 < £ < oo, with equal density of particles, p\ = pi = P- While transitions 
between a uniform (U) phase and a strip (S) phase are expected as p is varied, a novel 
droplet phase (D) which is not present in previous models |7|, [| apppears in between 
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(see Fig. 1). The three phases are characterized by different symmetries, with the 
localized droplet breaking translational invariance along both spatial directions and 
drifting steadily forward along its symmetry axis. That it is not an unstable or 
metastable state is revealed by the behavior of an order parameter, (ft (essentially the 
Fourier amplitude of the density profile at a smallest wavevector[||]), which varies 
from one (compact strip) to zero (perfectly uniform). Fig. 1(a) depicts 4> versus 
time (t) in the steady state for L = 48, £ = oo, showing frequent jumps among 
the three phases. In fact, a histogram analysis shows that with those parameters 
the system is at a three-phase coexistence, i.e., at a triple 'point. Deviating from 
this point, we find: (a) At slightly different density, the U or S phase becomes 
dominant, (b) For smaller bias £ with fixed L, <f>(t) plots and histograms show that 
the triple point splits into two transitions, e.g., S^D at psd, then D— >\J at pou-, a s 
p decreases, with gradually widening gap (psd — Pdu = — 0.02 as £ = oo — 0.5, for 
L = 48). Diminishing first order characteristics — hysterisis and finite jumps in eft 
and currents — suggest that the transitions eventually turn continuous, consistent 
with the observation that psd grows with decreasing £, making the droplet longer 
and fuzzier until no longer distinquishable from a strip, (c) For larger L at fixed £ , 
hysterisis becomes more pronounced while the transition points shift systematically 
downward. 

These results remind us of a pure substance which has a triple point and a first 
order phase boundary ending at a critical point Jl^]. Another possibility is to have a 
first and a continuous order phase boundary joining at a tri-critical point. Despite 
the obvious interest, however, we have not been able to map out the phase diagram 
completely. One difficulty is due to the above-mentioned decreasing transition den- 
sities for increasing system sizes, casting doubts on the survival of the transitions 
as L — > oo. At the same time, growing hysterisis makes it harder to locate the 
transitions precisely. To understand these crucial finite-size effects, we now focus 
our attention on a continuum theory. 

Mean-field theory: We begin by defining P^ n f n it...(x,y,t) as an equal-time, joint 
probability at time t for the locate state at sites (x, y), (x, y + 1), (x, y + 2), • • • to 
be n, n' , n", ■ ■ ■ where n = 0, 1 or 2 as above (-P^ n / n »... likewise). Since particles only 
hop to nearest neighbors, it is easy to see, for P n (x, y, t + 1) — P n (x, y, t) = P n , that 

\p\ = P [P?o(x,y-i)-PUx,y)] + <iiPoi(x,y)-Poi&y-i)} 

+ r[Pf (x - l,y) - Pf (x,y) + P^(x,y) - P \(x - l,y)}. (1) 

The 1/2 ensures the same time scale as in simulations, and lattice constants 1. 
P2 obeys a similar equation. An attempt to close this set of equations generates a 
hierarchy in the usual way. To proceed, we instead adopt a mean-field approximation 
in the following sense: 

PL? fav) « P n (x,v)P n ,(x,y + l)[l + A y nn ,(x,y)}, (2) 

with the asymmetric local correlations induced by biases accounted for by A nn /'s 
(similarly for P^ n i). Clearly A^ 2 > 0, A| x < etc. Only those with n / n' are kept, 
as A nn <C 1 is expected and well confirmed by simulation. 

By taking the (naive) continuum limit of Eq. ([!]) via Taylor expansions, assuming 
smooth functions, we obtain from Eq. (|l|) the continuity equation p n = — V • J n for 



3 



density p n , with the current 

Ji = xE lx pip 2 + y{Ep Q pi - Ei y p 1 p 2 ) 

-n \ ®Pi dp 2 , u .dpi dp 2] , , 

~ XD±[{1 ~ P ^ + P1 ^ ] ~ yD * [(1 ~ P2) ^ + * fl^' (3) 

where the parameters are completely specified by the microscopies: E = 2{p — q), 
E lx = 2r(Afj - Af 2 ), E ly = 2(pA y 12 - qA y 21 ), D\\=p + q, and D ± = 2r. J 2 can be 
obtained from Ji by switching the labels 1 2, and x *-* y. Before we continue, note 
that the terms proportional to E's and D's represent physically drifts and diffusions 
respectively. There is an induced drift normal to the bias in spite of symmetric 
transverse jump rates. Similar but much weaker corrections in the diffusive pieces 
due to A's are dropped in @. In the case of p = q, or p\ = 0, or p 2 = 0, the 
stationary state is a uniform product measure and Eq. (|2|) becomes exact with A = 0. 
Consequently, it is necessary to have two species driven to produce all the nontrivial 
behavior. Contrary to similar equations derived from entropy consideration 0, the 
microscopic origin of each term is made explicit. 

Recall that there are three different steady states in simulations. While the 
uniform solution is trivial to find, the droplet requires solving the full 2D problem. 
Since the droplet phase occupies an extremely narrow range in p (i.e., psd ~ Pdu), 
for practical purposes we may focus on the uniform and the strip phase. By setting 
J n =const, it is straightforward to show that the strip profiles satisfy 
1 l 

—Pi = Pi(l~pi+P2) (eipip 2 + C) (4) 

£o Po 
1 i 

Pi = P2O- + P1-P2) (eiPiP2 + C) (5) 

£0 Po 

for equal p\ = p 2 = p, where £0 = y/2E/(Dn + D±) = 2\/2(p — q)/(2r + p + q), 
ei = (E lx + E X y)/E = [pAf 2 + rAfi - qA y 21 - rAf 2 ]/(p -q),f/ n = d Pn /du, with the 
u axis perpendicular to the strip, at a polar angle of 135°. C = \[2J jE is a reduced 
current along +u, where J = Ji • (y — x)/y2. u and J have the same units as in 
simulations. 

Using the data for A's, we find a remarkable result: e\ ~ 1 — 2p = po for arbitrary 
£, L and p in the uniform phase (i.e., all data fall on one universal line). While 
apparently A — ► (so does £1) as po — > 0, we have yet an explanation for the full, 
universal behavior. In line with a mean-field approach, we assume it remains to 
hold locally and replace £1 by po(u) in Eqs. (Q) and (|5|). Agreement between the 



predictions and simulations provides a posteriori justification! 11]. Then Eqs 



and (|) contain no adjustable parameter and imply the scaling form 

p n (u,J,L,£) = f n (ue ,J/E,Le ). (6) 

Unlikely to be obtained in closed forms, the functions / n 's are determined by nu- 
merically solving the coupled Eqs. (Q) and (|5|). One stringent test of the theory is to 
compare these functions with those in simulations. Fig. 2 shows a typical test of not 
only f n but also the scaling in L. With no adjustable parameter, the agreements 
are impressive verifications of our approach. Deviation becomes noticeable only in 
the small Leo regime. Furthermore, by (y/2/L) J dup n (u) = p(J'/E,Leo), or 
equivalently 

C = V2J/E = g(p,Le ). (7) 
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This implies the entire current-density relation (along with the transitions) is the 
same for different L and microscopic £ , provided Leo is fixed. Fig. 3 confirms 
this prediction. The p dependence of g also compares favorably with simulations, 
except at small Leq. This is consistent with the observed transitions gradually 
turning continuous at small biases, as fluctuations ignored in the equations become 
significant. 

Eq. (0) allows an investigation of the thermodynamic limit. The transition point 
p c (L, £) is taken to be the minimum p allowed for the strip solution of a given Leq 
(i.e., at dC '/dp = oo; cf. Fig. 3). It is a function of Leo alone, hin^^oo p c (Leo) then 
decides whether the transition survives or not. Fig. 4 shows, for both simulations |ll 
and theory, that p c ~ (Leo) -0 ' 82 ^; and strongly suggests p c — > as Leo °o. 

These results have far reaching implications. If lim^Q^oo p c = indeed holds, 
then finite-density transition will occur generically only for infinitesimal bias (i.e., 
£ cx 1/L). It is hasty, however, to conclude that the transition would actually be 
extinct for finite £. Seldom found in equilibrium, this situation is by no means unfa- 
miliar in non-equilibrium steady states. A case of intense interest is the well studied 
ID limited local sandpile models of self-organized criticalityjl^]. Criticality of such 
models is controlled by the trough which has density pt vanishing as L -1 / 3 [14], but 



the number of troughs Lpt —* oo as L — > oo. Analogously, our model has the strip 
width £ ~ Lp c — > oo but £/L — > (hence £ is a mesoscopic length) as L — > oo, so 
that, like the sandpiles, the transition remains well-defined for all We suspect 

other models have similar features, but the corresponding analysis is hindered by 
anisotropics^ or by the lack of ergodicity||. 

To conclude, we have studied systems governed by asymmetric diffusion and two- 
species exclusion as simple models exhibiting non-equilibrium phase transitions. Our 
primary results are: (a) There exists multiple phases in spite of their apprarent sim- 
plicity; (b) a quantitatively successful theory with known microscopic dependence 
has been developed, which could be improved to handle fluctuations by adding noise 
terms, and extended to more complicated models or related problems, such as with 
inter-particle interactions; and (c) the survival of the associated transitions in the 
thermodynamic limit is addressed for the first time, and found to be parallel to 
certain class of model sandpiles. Finally, open questions such as the mechanism of 
stability for the droplet, details of the phase diagram and issues of critical exponents 
at small biases are of significant interests to the general study of non-equilibrium 
phase transitions. 
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Note added: At the completion of the manuscript, we received a preprint from 
Vilfan et al. [16], who obtained mean- field exact solutions for the strip phase in the 



oppositely biased two-species model [0]. Their results agree qualitatively with ours, 
including a tentative tri-critical point, but the large-L limits have not been explored. 
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Figure Captions 

1. (a) Order parameter at three-phase coexistence for L = 48, £ = oo, 
p = 0.1623: large, intermediate and small value corresponds to strip, 
droplet and uniform phase, respectively; (b) a typical strip for L = 128, 
p = 0.076; and (c) a typical droplet for L = 128, p = 0.074. Light (dark) 
particles drift to the right (top). 

2. Finite-size scaling of density profiles at fixed Lsq = 23.65 (or fixed C) 
and p = 0.36: L = 32(D), 48(A), 64(o), 128(x), 256(dashed line), and 
theory (solid line) [cf. Eq. (§)] with no free parameter. 

3. Scaling of current [Eq. (|7|)]. Thin line is p(l — p)(l — 2p) for uniform 
phase [from (|J)]. Branches are for strip phase at various Leo'- (a) 6.445, 
L = 48(o); (b) 14.81, L = 32(D), 48(A), 64(y); (c) 23.65, L = 48(«), 
64(A), 128(D), 256(x); (d) 33.84, L = 48(y); (e) 45.26, L = 48(D); (f) 
75.43, L = 80(y)- Note confirmed scaling in (b) and (c). Heavy lines 
are theoretical predictions for (b) to (f). 

4. Large-L behavior of transition density: theory [cf. Fig. 3] (A); data(» 
with error bars); linear fit (dotted line). 
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